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Quantum states can be stabilized in the presence of intrinsic and environmental losses by ei¬ 
ther applying active feedback conditioned on an ancillary system or through reservoir engineering. 

Reservoir engineering maintains a desired quantnm state through a combination of drives and de¬ 
signed entropy evacuation. We propose and implement a quantum reservoir engineering protocol 
that stabilizes Fock states in a microwave cavity. This protocol is realized with a circuit quantum 
electrodynamics platform where a Josephson junction provides direct, nonlinear coupling between 
two superconducting waveguide cavities. The nonlinear coupling results in a single photon resolved 
cross-Kerr effect between the two cavities enabling a photon number dependent coupling to a lossy 
environment. The quantum state of the microwave cavity is discussed in terms of a net polarization 
and is analyzed by a measurement of its steady state Wigner function. 


An unavoidable adversary in quantum information sci¬ 
ence is decoherence. A large scale quantum computer 
must implement error correction protocols to protect 
quantum states from decoherence [1]. A first step toward 
fault tolerant quantum error correction is the stabiliza¬ 
tion of a particular quantum state in the presence of deco¬ 
herence [2]. One such implementation uses gate based ar¬ 
chitectures with measurement and feedback [3-9] for the 
correction of quantum errors. An alternative approach to 
active quantum systems is quantum-reservoir engineer¬ 
ing (QRE) [10-13] which harnesses persistent, intentional 
coupling to the environment as a resource. Both cases 
require entropy removal yet only QRE employs environ¬ 
mental losses as a crucial part of their protocols. QRE 
does not require an external feedback with calculation 
since the Hamiltonian interactions are designed a priori 
to determine the final state avoiding uncertainty induced 
by the quantum-classical interface. In addition, QRE is 
less susceptible to experimental noise [14] and in some 
cases thrives in a noisy environment [15]. 

QRE has been demonstrated in macroscopic atomic 
ensembles [16], trapped atomic systems [17], and su¬ 
perconducting circuits [18-20]. Circuit quantum elec¬ 
trodynamics (cQED) systems are an attractive platform 
for QRE due to the experimental freedom to design 
strong interactions between superconducting qubits and 
microwave cavities [21]. Interactions between a super¬ 
conducting transmon qubit and a microwave cavity have 
demonstrated qubit-photon entanglement [21], creation 
of quantum oscillator states [22, 23], and quantum non¬ 
demolition measurements of an oscillator [24]. Investiga¬ 
tions using three dimensional waveguide cavities resulted 
in increased coherence times [25-27] allowing the obser¬ 
vation of novel quantum phenomena such as the single 
photon self-Kerr effect [28], a protocol that creates arbi¬ 


trarily large Schrodinger cat states [29], realtime parity 
monitoring of the decay of cat states [30], and a protocol 
that confines the state of a cavity to a quantum manifold 
[31], yet no demonstration of a nontrivial cavity state 
QRE protocol exists. 

In this Letter, we demonstrate the single photon re¬ 
solved cross-Kerr effect between two superconducting mi¬ 
crowave cavities which is a new regime of cQED. This 
nonlinear coupling causes an excitation in one cavity to 
change the resonance frequency of the other cavity by 
more than their combined linewidth. While the state 
dependent shift between a qubit and a cavity has been 
previously observed [32] in this Letter we present the 
observation of the state dependent shift between two mi¬ 
crowave cavities. In this work a transmon is used to 
introduce nonlinearities to the cavities and for tomogra¬ 
phy. This new regime of cQED is used to realize the first 
cQED QRE protocol that stabilizes quantum states of a 
microwave cavity. Here we demonstrate a protocol that 
stabilizes a one photon Fock state of a microwave cavity. 
Since the storage cavity is restricted to its first two en¬ 
ergy levels, this process can be described as a population 
inversion and as an effective negative temperature. This 
protocol could be extended to higher photon states of 
the microwave cavity by including more CW drives. The 
single photon resolved cross-Kerr is necessary for a QRE 
protocol that stabilizes cat states of an oscillator [33] and 
may be used as a cavity-cavity entangling operation. 

Within a cQED framework, we model our system as 
two harmonic oscillators coupled to a nonlinear oscilla¬ 
tor [34]. The most nonlinear oscillator in our circuit is the 
transmon whose nonlinearity originates from the Joseph¬ 
son junction with an inductance that is nonlinear with 
respect to the flux across it. This system is well described 
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by the following Hamiltonian [25, 29, 35] 

+ tUcC^c 

_ 

2 2 2 
- Xgsa'^'ab^b - x^^a'l'ac'l'c - Xscb'^bc^c (1) 

The subscripts used in the Hamiltonian are ‘g’ for the 
transmon qubit, ‘s’ for the storage cavity, and ‘c’ for 
the cooling cavity. On the first line are dressed angu¬ 
lar frequencies denoted by uji. The second line contains 
self interaction Kerr type terms, called anharmonicities, 
of the modes denoted by Ai. On the final line are the 
state dependent shifts, between modes. Since the 
state dependent shift, to fourth order, is proportional to 
the geometric mean of the anharmonicities of the modes, 
Xij = 2^yAiAj [35], a strong, dispersive interaction be¬ 
tween modes requires an appreciable anharmonicity for 
each mode. 

The experimental setup consists of two microwave cav¬ 
ities [28, 29] machined from high purity (99.99%) etched 
aluminum [26]. Between the two cavities is a trough 
where the transmon is positioned to couple the cavities. 
Couplers are used to set the dissipation rate of the cool¬ 
ing cavity Kc, while the coupler to the storage cavity is 
weakly coupled allowing direct excitation of the storage 
cavity. The transmon is fabricated [36] on high purity 
c-plane sapphire and is coupled to each cavity by an an¬ 
tenna which enables a strong, dispersive interaction be¬ 
tween the qubit and each cavity [28] Xqs/27r(21.1 MHz) > 
'=s/27r(65 kHz) and Xgc/27r(4.9 MHz) > «c/27i-(1.7 MHz). 

To measure the storage cavity anharmonicity, we use a 
single CW drive to perform spectroscopy. Fig. 1(a). Us¬ 
ing a large amplitude drive which power broadens the 
/s,o_>.i {fi = ‘^i/27r) transition, we observe the two pho¬ 
ton transition with frequency . The detuning cor¬ 

responds to half the anharmonicity, Hg, of the storage 
cavity and we infer an inherited cavity anharmonicity 
^s/27r = 4.0 MHz. Following the same method, we de¬ 
termine that the cooling cavity has an anharmonicity 
^c/27r = 300 kHz. To measure the state dependent shift 
between the two cavities we first perform a 5 ns displace¬ 
ment pulse on the storage cavity, then spectroscopy on 
the cooling cavity and finally high-power readout [37]. 
Shown in Fig. 1(b) is a spectroscopy measurement of the 
cooling cavity for a displacement n « 1.5 of the storage 
cavity. Discrete spectral peaks for up to three photons in 
the storage cavity are visible. From this we infer a state 
dependent shift = 2.59 ± .06 MHz and observe the 
first single photon resolved cavity-cavity cross-Kerr. 

The measured Hamiltonian parameters lend them¬ 
selves well to a cQED QRE protocol that stabilizes Fock 
states. The first requirement for this protocol is that 
the cavity in which Fock states will be stabilized is more 
anharmonic than its natural linewidth, A^ > Ks so that 
individual transitions may be selectively driven [Fig. 2(a) 
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FIG. 1. Storage and cooling cavity spectra, (a) Spectroscopy 
is performed on the storage cavity with a single CW drive. 
With a large amplitude drive, we observe the two photon tran¬ 
sition, . From this measurement we infer the location of 

the /s,i -»2 transition (red dashed line) and determine its de¬ 
tuning from the /s.o^i transition as 4.0 MHz which we define 
as the anharmonicity of the storage cavity, (b) A 5 ns dis¬ 
placement pulse, whose amplitude gives n « 1.5 in the stor¬ 
age cavity, enables the observation of a single photon resolved 
cross-Kerr between the two cavities, = 2.59 ± .06 MHz. 


left]. A second requirement is a state dependent shift be¬ 
tween the two cavities that is larger than both of their 
linewidths, Xsc > [Fig. 2(a) right]. In Fig. 1 we 

see that these requirements are met. However, this pro¬ 
tocol is most successful when the lifetimes of the storage 
cavity and the cooling cavity are quite different Kc'A> Kg. 
In Fig. 1(b) the decay rates of the cavities are compara¬ 
ble, ('^c/kb « 4). We alter the ratio of lifetimes, while 
maintaining the necessary requirements Ag > Kg and 
Xsc > ^ Kg, between the two cavities by increas¬ 

ing the coupling strength to the external environment of 
the cooling cavity resulting in a new ratio of 25. 

Shown in Fig. 2(b) is a QRE protocol that stabilizes 
a one photon Fock state in the storage cavity. This pro¬ 
tocol is conceptually similar to the protocol used in Ref. 
[19] which stabilized the ground state of a qubit tensor 
product with a coherent state of a cavity. Although we 
stabilize the ground state of the storage cavity, we also 
use this protocol to stabilize a one photon Fock state. 
Due to the anharmonicity of the storage cavity, a CW 
drive, Us, can be applied to the /s,o->.i transition. This 
drive is an induced Rabi rate between vacuum and a one 
photon Fock state. Concurrently, with Us, a drive with 
strength Uc is applied detuned by a cross-Kerr from the 
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cooling cavity. This drive is resonant provided that there 
is exactly one photon in the storage cavity. Once reso¬ 
nant, the conditional drive displaces the cooling cavity 
to a coherent state determined by the amplitude of the 
drive. When a photon decays from the storage cavity 
Oc, is no longer resonant and the cooling cavity quickly 
decays to vacuum. Once back to the ground state, the 
storage cavity is resonant with the drive fig. This pro¬ 
tocol reaches its steady state solution in a time governed 
by the decay rate of the cooling cavity. The steady state 
population in the one photon Fock state of the storage 
cavity will be determined by its decay rate, Ks, and the 
stabilization rate, The stabilization rate is defined 
as the rate at which the system is returned to the target 
state when a photon decays from the storage cavity. Us¬ 
ing a simple four state model we expect that to achieve a 
99% one photon Fock state in the storage cavity, a mini¬ 
mum ratio of lifetimes between the two cavities of 300 is 
required see supplementary material. 

This protocol requires both the frequencies of the two 
microwave drives and their amplitudes be chosen appro¬ 
priately. From a full simulation of the Linblad master 
equation as well as our experimental observations, we 
hnd optimal performance when Us ~ ^c- We determine 
the drive power applied to the cooling cavity through a 
power dependent dephasing measurement of the trans- 
mon qubit applied roughly at one cross-Kerr detuned 
from the cooling cavity [38]. 

The experimental implementation begins with CW 
drives applied simultaneously to the storage and cool¬ 
ing cavity for a duration of 200 k“^ which is twenty times 
longer than the time necessary to reach steady state [19]. 
To measure the photon population in the storage cav¬ 
ity, we stop the drives, wait for photons to decay from 
the cooling cavity, and apply conditional qubit tt pulses 
to determine the photon number in the storage cavity 
[24, 28]. We vary both the drive strength and frequency 
applied to the cooling cavity while maintaining optimal 
parameters for the storage cavity. 

We plot the steady state polarization, p = , of 

the storage cavity after running the protocol in Fig. 3(b). 
P{n) corresponds to the probability of having n pho¬ 
tons in the storage cavity. Due to the selectivity of the 
drive, Ug, the storage cavity is limited to its first two 
Fock states. We confirm this by measuring populations 
for the two and three photon Fock states. When Uq 
is driven at the zero photon peak of the cooling cavity 
we observe p = 0.95 demonstrating that storage cavity 
is overwhelming in the zero photon Fock state despite 
the induced Rabi drive on the storage cavity. However, 
as the drive power and frequency applied to the cooling 
cavity are varied, steady state stabilization of a polariza¬ 
tion inversion occurs corresponding to a predominantly 
one photon Fock state in the storage cavity. This pop¬ 
ulation inversion is a purely quantum effect and can be 
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FIG. 2. Ideal cavity spectrum and Fock state stabilization 
protocol, (a) Left: sketch of idealized storage cavity spec¬ 
trum. The storage cavity must have unequal energy levels 
spacing (hAs), inherited from the coupled qubit, to selectively 
drive storage cavity transitions. On the right is the idealized 
cooling cavity spectrum. The frequency shift of the cooling 
cavity due to photons in the storage cavity, the cross-Kerr 
(Xsc), must be larger than either cavity linewidth to selec¬ 
tively drive this transition, (b) Energy level diagram for the 
coupled cavity-cavity system tracing over the qubit state. As¬ 
cending vertically are excitations in the storage cavity while 
to the right is increasing number of excitations in the cool¬ 
ing cavity. A microwave drive, Us, is applied on the storage 
cavity so that population only oscillates between vacuum and 
the first Fock state of the storage cavity. Simultaneously a 
drive, Dc, is applied on the cooling cavity such that it is 
resonant provided there is exactly one excitation in the stor¬ 
age cavity. Once resonant, the cooling cavity is pumped to a 
mean photon number set by the strength of the drive. The 
autonomous loop of this protocol is closed by cavity decays, 
decaying arrows, returning the population to |0,0) allowing 
the preparation to be repeated. 


described as an effective negative temperature according 
to [39]: 

2kB tanh'^ (p) 

Where h is Planck’s constant and fee is Boltzmann’s con¬ 
stant. From this we infer that our steady state solu¬ 
tion corresponds to an effective negative temperature of 
—0.77 ± 0.06 K in equilibrium with the storage cavity. 

In Fig. 3(d) we see the time dynamics of this proto¬ 
col where the initial polarization is unity then chang¬ 
ing as a function of time to its steady state value of 
p = —0.26 ± 0.04. Plotted on top of the data is a full 
simulation of our driven dissipative system where we find 
excellent agreement in our time dynamics [40]. From 






















4 


detect N in storage 



Stabilization time (Ts) ^ ^ 



Wigner Tomography 




-1.5dB 
-8.5 dB 
Theory 


0 12 3 

Normalized detuning (A/zg^) 


c 

g 

_ro 

o 

Q. 


1.0 


—^-p^l 1 1 ' 1 M 

0.0 


region of T<0 ^ 

1.0 


-^-Lr 1 1 1 1 1 1 1 id 


0.5 1.0 40 80 120 160 200 

Stabilization time (ps) 


FIG. 3. Storage cavity polarization, (a) The Fock state sta¬ 
bilization protocol described in Fig. 2(b) is applied for a du¬ 
ration, Ts, followed by a 300 ns wait to evacuate excitations 
from the cooling cavity, a photon selective tt pulse is then 
performed on the qubit determining the probability of each 
photon state of the storage cavity up to three photons, (b) 
Storage cavity state polarization as a function of drive ampli¬ 
tude and frequency. The frequency of the cooling cavity drive 
is plotted as A = /q — /d, and normalized by the cross-Kerr, 
Xsc- As the frequency of the drive applied to the cooling cav¬ 
ity is brought in resonance with the first photon peak of the 
storage cavity ^ « 1 the protocol stabilizes the first Fock 
state of the storage cavity. The inset is a simulation plot with 
the same axis and color scale as the experimental result, (c) 
Linecuts for a weak drive power and a drive power resulting in 
a polarization inversion, (d) As the duration of the stabiliza¬ 
tion protocol is varied the polarization of the storage cavity 
alters and for infinite time reaches its steady state solution. 


the four state model, we would expect a polarization of 
p = —0.47. This value is within a factor of two of both 
what is measured experimentally and extracted from a 
full simulation of the Linblad equation. Through simula¬ 
tion of the full Linblad master equation we find that the 
limitation in polarization inversion is due to the finite 
ratio of lifetimes. Physically realized two cavity systems 
exist with a ratio of lifetimes being a factor of one thou¬ 
sand [301. With that ratio we expect that in steady state 
P(l) > 0.99. 

Although much of the discussion has framed the stor¬ 
age cavity in the language of spin systems it is still an 
oscillator. To demonstrate the oscillator nature of the 
storage cavity in Fig. 4 we perform cavity tomography 
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FIG. 4. Wigner tomography of stabilized steady state of the 
storage cavity, (a) The previously described stabilization pro¬ 
tocol is used to reach the desired steady state. Then Wigner 
tomography is performed on the state of the storage cavity, 
(b) Left: Measured Wigner function for the steady state of 
the storage cavity which is a statistical mixture of an A = 1 
and N — 0 Fock state. Right: Simulated steady state of the 
protocol, (c) Linecuts along Im(Q:) and Re(a) for the mea¬ 
sured Wigner function and the simulated steady state Wigner 
function. Although not a pure N — 1 Fock state of the stor¬ 
age cavity our long term solution does have negativity in the 
Wigner function indicative of a quantum state. 


measuring generalized Husimi Q functions, Qn{o) = 
■TT~^\{N\D-a [28], up to A = 3 Fock state of the 

storage cavity, D-a is the displacement operator, and lb 
is the final state and infer the Wigner function by adding 
and subtracting the even and odd measured Q functions. 
We compare these results to the Wigner function of a 
simulation of the steady state solution to the protocol. 
Our results are explained in terms of a harmonic oscil¬ 
lator picture with the steady state of the storage cav¬ 
ity in a statistical mixture of P{0) = 0.37 ± 0.03 and 
P(l) = 0.63 ± 0.02. In Fig. 4 c statistically significant 
negativity in the Wigner function is observed. 

In conclusion, we present the first single photon re¬ 
solved cross-Kerr effect between two cavities. We used 
this system to implement a cQED QRE protocol that 
stabilizes Fock states in a superconducting microwave 
cavity. We demonstrate one such instance, stabilizing 
an A = 1 Fock state, by presenting its reconstructed 
Wigner function to demonstrate the oscillator nature of 
the storage cavity. This protocol can be extended to 
higher photon numbers of the storage cavity by including 
more selective microwave drives at the different transi- 
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tions of the storage cavity. Our steady state polarization 
inversion corresponds to p = —0.26 ± 0.04 which we map 
to the storage cavity being in equilibrium with a bath 
of T = —0.77 ± 0.06 K. Our protocol is limited by in¬ 
duced spontaneous emission to the environment. Future 
implementations would benefit from a Purcell filter and 
increased nonlinearity in the cQED system. 
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I. SYSTEM AND HAMILTONIAN PARAMETERS 


The sample is a vertical transmon mounted in a high purity, etched aluminum cavity. This sample is secured in a 
dilution refrigerator whose base temperature is 20 mK. The filtering, magnetic protection and mounting are identical 
to[l]. Readout is performed via high power readout of our sample using the high power peak of the cooling cavity [2]. 
Our full Hamiltonian to fourth order in the junction flux is: 



Storage 

cavity 


transmon 


cooling 


FIG. 1. Picture of sample 
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Term 

Value 

Units 

^9/277 

7249 ± .013 

MHz 


8493 ± .017 

MHz 

‘^c/ 27 r 

9320 ± .024 

MHz 

^9/277 

26 ± .3 

MHz 


4.0 ± .1 

MHz 

"^=/ 27 r 

300 ± 79 

kHz 

Xqs/2-n: 

21.1 ± .05 

MHz 

Xq <^ j 2 TC 

4.9 ± .04 

MHz 

XstLj2-w 

2.59 ± .06 

MHz 


In this section we also include the decay rates of the cavities for the measurements of figure 2 of the main text 
where Ks = 7.5 kHz and Kc = 120 kHz. The measurements done in figure 4 of the main text for the stabilization 
protocol are Ks = 65 kHz and Kc = 1.7 MHz. 


II. FOUR LEVEL MODEL 

To develop intuition as to how well the Fock state stabilization protocol should work we will look at the simplified 
case of a four state system. Starting with the state in bottom left as state A we will label states clockwise as B, C, 
and D. Our intention is to stabilize state C. To determine how well we can stabilize this state we must determine the 
ratio between the rate that state C decays to state D which we will call The other rate that we must determine 
is the rate at which our system corrects the C to D decay process and returns the system to state C. We will call this 
the stabilization rate and refer to it as k^. This rate will be a combination of the decay rate from D to A, kda, the 
rate to be driven from state A to state B, ^Iab, and the rate at which state B is driven to state C, kbc- Since states 
C and D serve as a proxy of an oscillator we will have the rate at which it decays, kda, being equal to the rate at 
which it rings up, kbc and just call this rate k. To recap: 


kda = Kbc = 


^BC 
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FIG. 2. Four state system with decay rates 


Now we must determine how k and f^AB are used to determine Kf. We identify that if either k or Qab are zero 
then K-j- should also be zero regardless of the value of the other rate. For this reason we add the rates in inverse: 
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The optimal choice is for 0.ab ~ ^ found through an eigenvalue calculation of the matrix modeling for the four state 
model. With this simplification we find: 


K 

«T=3 

In the absence of decay from the target state we would expect unit hdelity in creating it. When accounting for decay 
we expect that the ratio of and K-f describes how far from unity the stabilized target state will be. For instance if 
the two rates are equal we expect that our stabilized state will be produced with a fidelity of 0.5. In equilibrium we 
have: 


K^Pil) = n^P{0) 


Solving for P(l) using that P(l) + P(0) = I we get: 


P(I) 


1 


1 + IT 


From this we see that if we would like to have a fidelity of 0.99 to the target state the would require a ratio between 
the decay rate of the cooling cavity, k, and the decay rate of the target state, of 300. From simulation we find that 
for our actual system that |. 


III. LINBLAD MASTER EQUATION AND SIMULATION 

We begin with just cavity terms in our Hamiltonian: 

^/h = to'sb'^b + to’cc'^c-^b'^^b^-— Xscb^bc'^c (2) 

We add in drives of the form Qs(b^ + b) for the storage cavity drive and Qc(c^ + c) for the cooling cavity drive. 
Where Qs and flc are complex drive amplitudes. By including these drives and entering the rotating frame of both 
cavities Ug = and Ug = result in detunings of Ag = ojg — for the storage cavity and 

Ac = Wc — idjfic for the cooling cavity. Our driven Hamiltonian is now: 

^/h = Agb'I'b + AcC^c - ^b^^b^ - - Xscb'^’bc'l'c + ris(b''’ + b) + Oc(c'^ + c) (3) 

With this driven Hamiltonian we can write the master equation for our system as: 

p = -*[H,p] + nsD\b]p+ KcD[c]p 

The above master equation is used in QuTIP 2.2 [3] with the steady state solver and the time dependent solver to 
produce the simulation results that appear in Figure 3 and Figure 4 of the paper. 
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